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ABSTRACT 


A procedure was developed to aid in the testing of a 
two-sided composite null hypothesis about the mean of a 
normally distributed random variable for situations where 
either the population variance is known or unknown. The 
procedure was designed to eliminate the requirement for 
iterative type solution techniques normally used in deter- 
mining the acceptance or rejection region of the subject 
hypothesis. This thesis provides guidelines, curves, and 
tables which will aid in testing a two-sided composite null 
hypothesis. Provisions were also incorporated into the 
procedure to permit testing of hypotheses about the relative 


displacement of the coefficient of variation from zero. 





TABLE OF CONTENTS 


4E INTRODUCTION ------------ --- ---- -- -- -- -- - - - - - -- -- - - 9 
A. BACKGROUND AND PURPOSE ------------------------ 9 
В. SCOPE ---------- -- -- ------ - - -- --- - - - - -- - - - - - - -- 9 
C. SYMBOL AND HYPOTHESIS TESTING CONVENTIONS 
PECULIAR TO THIS THESIS ----------------------- 1 
Il. RESULTS AND DISCUSSION ---------------- - - -- --- - -- -- 12 
A. CASE OF KNOWN POPULATION VARIANCE ------------- 12 
1. General Remarks --------------------------- 12 
2. Procedure Description --------------------- 13 
a. For Controlling a Error Level --------- 13 
b. For Controlling 8 Error Level --------- 14 
c. For Estimating Tail Probability E -- 14 
3. Range of Parameter Values Considered ------ 17 
4. Example Applications ---------------------- 18 
B. CASE OF UNKNOWN POPULATION VARIANCE ----------- 20 
l. General Remarks --------------------------- 20 
2. Procedure Description --------------------- 21 
a. For Testing That |u- u < rg ---------- 22 
b. For Testing that CV < Ne E 22 
c. For Estimating Tail Probability (P4) -- 23 
3. Range of Parameter Values Considered ------ 2 
4. Example Applications ---------------------- 28 
APPENDIX A: CURVES FOR DETERMINING a OR 8 CRITICAL 
VALUES AND TAIL PROBABILITY (Pm) ----------- 30 





APPENDIX B: TABLES OF a OR B8 CRITICAL VALUES AND B, ------ 47 


T 
APPENDIX C: PROCEDURE DERIVATION ------------------------- 94 
APPENDIX D: SOLUTION TECHNIQUES AND COMPUTER PROGRAMS ----104 
APPENDIX E: LIST OF REFERENCES ---------------- ----------- 116 
INITIAL DISTRIBUTION LIST --------------------------------- IT 





E C 


L'"R 


να | 


С* 


EV 


CV 
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SYMBOLOGY 


Probability of a Type I error. 

Probability of a Type II error. 

Critical value defining the rejection or acceptance 
region centered around Ho (Oo known case) or zero 


(o unknown case). 


The lower bound (Cr) and upper bound (Cp) on the 
critical region. 


| 
х. 
del 


n 





The true population expected value (mean). 


A specified value of the true population mean. 


. Specified lower (d. ) and upper (dp ) limitsrfor 


the true population mean (L). 


Population variance (o = population standard 
deviation). 
Sample variance used to estimate ο. specifically; 


(x. - ays 


sample size. 


С Уп 


σ 


ὃ yn 


σ 
Coefficient of variation К 
ТТ) 
A specified value of the CV. 


The null hypothesis of interest. 





DF 


Ф [х] 


The alternative hypothesis to Ho 
The tail probability that a value equal to or 


more extreme than the observed X would have 
resulted if Ho were in fact true. 


Desired maximum allowable ο of the true 
mean from some specified value (и in terms of a 
portion of the population Ns deviation 
(i.e., |u- ug| < ro). 
The statistic 

n(x = u) 


Degrees of freedom. 


The non-centrality parameter of the non-central 
Todgrstributlone A =r n- 


Normal cumulative distribution function. 








I. INTRODUCTION 


A. BACKGROUND AND PURPOSE 

A procedure was developed to aid in testing a two-sided 
composite null hypothesis about the mean of a normally dis- 
tributed random variable. The intent was to eliminate the 
need for iterative solution techniques when solving for 
critical values manually. In situations where computer 
assistance is inappropriate or unavailable, there is a tendency 
for the analyst to shift from the interval test requirement 
to the computationally easier two-sided simple hypothesis 
(i.e., that the mean is equal to the midrange value of the 
interval). By employing the procedure proposed in this thesis, 
the time required to test a two-sided composite null hypothesis 


by hand can be significantly reduced. 


B. SCOPE 
Procedures were developed for the following test requirements. 


IU cr Testing Thdtedeses n scd 


L R 
For situations where the population variance is known, 
a method was developed for testing that the population mean 


JE 


is between two values (ч, and ав 


2. For Testing That ШЕШ < ro 
For situations where the population variance is unknown, 
a method was developed to test that the absolute value of the 


deviation of the population mean from some specified value (н) 





is less than or equal to a specified portion (г) of the 
population standard deviation. By a slight modification, 
this procedure can be used to test if the coefficient of 
variation is between zero and some specified value. 

Each procedure consists of using relevant situation 
parameters (i.e., б, п, 6, a, etc.) to determine the critical 
region using EE interpolation curves or tabulated data. 
The procedure was designed to be versatile with respect to 
using Type I (a) and Type II (8) error level control, or 
using tail probabilities (Р) for drawing statistical 
inference. 

To simplify use of the procedure, the paper is 
organized so that the main body contains only 1) minor 
theoretical context, 2) detailed procedural guidelines, and 
3) example applications. The detailed derivation of the pro- 
cedure is provided in Appendix C for readers so inclined. 
Appendices A and B provide curves and tables used to find the 
associated rejection or acceptance regions and the tail area 
probability. The curves are provided to simplify critical 
value interpolation for what is considered to be the majority 
of applicable parameter values. However, the total breadth 
of parameter ranges considered in this paper is covered in 
the tabulated data of Appendix B. It should be noted that 
all explanatory figures in the main body of this paper are 


Schematic in nature and do not represent scale drawings. 
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C. SYMBOL AND HYPOTHESIS TESTING CONVENTIONS PECULIAR 
TO THIS THESIS 


Whenever possible, symbol and concept conventions which 
conform to standard hypothesis testing terminology were 
employed. However, because of requirements which are peculiar 
to this procedure, it was convenient in some instances to 
use terminology that may be new to the reader. The following 
symbol and concept conventions apply throughout this Thesis: 

l. The Null Hypotheses (Ho) will always be bounded 
intervals such as [m m] <0 or | < ro. All alterna- 
tive hypotheses are the complements of such intervals. 

2. The level of significance is customarily a pre- 


assigned upper bound on the size of the critical region 


(i.e., that subset of the sample space reserved for rejection 
ОЁ Н). This will hereafter be referred to as a error level 
control. 


3. By reversing the role played by the null and alterna- 
tive hypotheses, it is possible to define an acceptance 
region for the originally formulated interval null hypothesis 
with a specified preassigned upper bound (8) on the size of 
a type II error. This procedure is defined as 8 error level 


control. 
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II. RESULTS AND DISCUSSION 


A. CASE OF KNOWN POPULATION VARIANCE 
l. General Remarks 
In order to test the null hypothesis that the popu- 
lation mean falls between two values (i.e., d, < u< дь), 


L 


an iterative solution of equation 1l is required. 


where, 


G is the known population standard deviation, 


C is the value of the statistic x which defines 
the critical region centered around 


NL. oe 
Ho 3 , 
a is the probability of a Type I error, 
ΝΞ 
ô is equal to που and 
Ф is the integral of the standard normal probability 


density. 

To determine a rejection region for the subject 
hypothesis, an initial value of C would be guessed and the 
right side of equation 1 calculated. A new C would then be 
used iteratively until the desired accuracy in the value of 
a had been obtained. This solution process is not difficult 
or time consuming with the assistance of a computer. However, 
when performed manually, such a process may require more effort 


than is warranted by the scope of the project. 
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The procedure developed in this thesis will provide 
the individual with a method for testing a two-sided com- 
posite null hypothesis which circumvents the previously des- 
cribed iterative solution technique. The procedure derivation 
is provided in Appendix C. In general, it consists of 


6/n 


normalizing equation l by reparameterizing the terms —— 


σ 
C/n T 


and — as the parameters "a" and C*, respectively.  Appro- 
priate values of C* are then obtained according to the pre- 
Scribed procedure which follows. 
2. Procedure Description 
The procedure for testing that dr SOS da is out- 
lined in Figure 1 (for a and 8 error control) and Figure 2 


(for estimating P This procedure was developed for 


T)* 
situations where symmetric critical regions are of interest. 
A summary of the procedure follows. 
a. For Controlling The a Error Level 

First calculate a - uE If the value of "a" 
is less than or equal to 1.7, the associated value of C* 
is obtained from Figure А-1 of Appendix A. The lower (C; ) 
and upper (Cp) bounds of the rejection region are then 


defined using equation 2. 


C*g C*g I R 
A Rog lm 








If the calculated value of "a" is greater than 
1.7, the relationship between "a" and C* becomes linear and 
equation 3 is used to define the rejection region (Κε 15 


obtained from Figure 1). If x < C, or x > C, then reject 
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) 


the hypothesis with a probability equal to gq of rejecting 
when in fact the hypothesis is true. 
(a + Kp)0 la + K.) o 


3 = ο ον ο ον = 
o) Sy Esc - ο ο ο E 


B: Por Controlling The B Error Level 

First calculate a = > If "a" is less than or 
equal to 3, the associated C* value is obtained from Figure 
Seem 2 5 a > 1-6) or Figure A-3 (1.6 < a < 3) of Appendix A. 
The lower (Cy) and upper (C4) bounds of the acceptance 
region are then defined using equation 2. 

If "a" is greater than 3, the relationship 
between a and C* becomes linear and equation 4 is used to 
define the $ error level acceptance region. If C, US Cor 
the hypothesis that the mean is between dr and da is accepted 
with a probability of falsely accepting equal to 8. 

el ج‎ К.) 0 (а - Κα) σ 
(4) C, 7 Mg 7 = and C, = u + чое 
с. For Estimating Tail Probability (Р) 

The exact value of the tail probability (Р) сап- 
not be calculated Since it depends on the actual value of the 
population mean (u). However, the procedure depicted in 


Figure 2 will provide an estimate of the maximum possible 


value of P_ for a corresponding value of x equal to or more 


T 
extreme than that observed. The procedure for estimating Р 
consists of first calculating the values a = S and 
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бит 


G + 


ar >" 5 


Go to Append.A Co to Append.A(Fig.2 
εία.]. TOT Ca l0 OF 19.3 
for 1.6«a«3). 





Calculate a= 


















Calculate " & C 
_ (a - Keo 
EOS /n 
( - Κε 
CR Uo" š Yn Е 


Calculate C at desired error 
level curve and intersection 
of vertical line projected 
upward from calculated value `a; 








> 















Calculate C, and Ca з 
* с“ 
! к σ 
Cp uot yn 





If x «C or Х>Ср then 


reject H (d, < u < dp) 


If Ск 
H (d < n) 
probability of falsely 
accepting 


u < d with an & 


with an a probability 
of falsely rejecting. 


APPROPRIATE Kr VALUES 


eve E 
d, á d.-Lower and upper test interval 
values. 


S(dp- d, ) 


-5(d, + dp) 


Defines lower bound of critical 
region. 


C,= Defines upper bound of critical 
region. 





Figure 1: Algorithm For Determining a or 8 Critical Regions 
for Known © 
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Yes 


B 2 
in P a 5 





ТЕ 


Go to Append. A Calculate PT as: Calculate Pr from: 
(Fig. 4 for 0<а<]1. 5 > ж-о, г. pt Zafra pr 
oe Fig. 5 for Pr l-da*C ]teia-C ], Рт= ФГа-С ] 


1.5<а<2.3) апа deter- where ὁ is the integral 
mine point where of the standard normal 
vertical line from density. 

'a' value and 
horizontal line 


through C" intersect. 
















Interpolate between 

mes to determine 
(Note: these curves 

E. plotted log 

and should be inter- 

polated accordingly). 


Figure 2: Algorithm For Determining The Tail Probability 
For o Known 
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Ск = 


(x - μο) ντι 
ف‎ These values are then used to estimate 


σ 


Р by employing the interpolation curves of Appendix A or 


by direct calculation using the equations provided in Figure 2. 
3. Range of Parameter Values Considered | 
The procedure itself is suitable for all possible 
parameter values. However, a restriction on the range of 
parameters covered by the curves and tabulated data was 
necessary. The range of parameter values was selected ina 
manner which would provide adequate coverage for the majority 
of conceivable real-world situations. 
a. Curves 
1. For a Error Level Control. Parameter ranges 
ο теш ате с.о, 025, 405, 075, .1,..15, .2, and .25) 
ana a(0 to 1.7). 
2. For B Error Level Control. Parameter ranges 
Ep cred'are B(.01, .025, .05, .-.075, .1, .15, .2, and .25) 
End a(.2to 3). 


3. For Estimating P Parameter ranges covered 


т” 
аге а(0 to 2.8) and C*(.04 to 5.2). 
b. Tables 
1. For a Error Level Control. Parameter ranges 
covered are %(.001, SOM ολο OS OSS CI CIS, a2, And 
#5) and a(0 to 3). 
| 2. For 8 Error Level Control. Parameter ranges 


ссе ае οοἳ οἱ, 025, —05,. 2075,91, .15,..2, and 


725) and a(0 to 5). 
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4. Example Applications 
a. Example 1 

In the production of a certain Structural element, 
a new size has been established as marketable. The expected 
size of the new element is to be between 16.5 and 17 cmin 
order to satisfy engineering and management requirements. 

For example, if the expected size is below 16.5 cm the 
rejection rate of unacceptable elements will be too high. 
Be the expected size exceeds 17 cm, more material is used 
than is necessary. 

For the first 25 units produced, the average 
size is 17.03 cm. It is known from previous experience that 
the element size is normally distributed with a standard 
deviation equal to .2. In addition, no functional relation- 
Ship between the mean and standard deviation of the element 
Size has been observed in the range of 10 to 20 cm. 

The decision maker is interested in determining 
if the hypothesis that the expected size is between 16.5 and 
17 cm can be rejected with a Type I error equal to .05. The 
test is NOU Cd with 6 equal to .25, Ho equal to 16.75, and 


251305) 
2 


a equal to or 6.25. Since a is greater than 1.7, 


the rejection region is defined using equation 3, specifically: 


EE (ло, 


L E = 16.434 


and 
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Cu zn ce 3158 5 17.066. 


Since x is neither greater than Cp or less than C the 


L’ 
hypothesis cannot be rejected with a Type I error equal to 
505. 

b. Example 2 

From a sample size of 81 shipping containers 

selected at random, it 1S necessary to know if the expected 
container weight (which appears to be normally distributed) 
1s between 75 and 75.5 kg. The sample average weight and 
variance are 75.22 and 1.44, respectively. A Type II error 
equal to .01 is desired. Since the sample size is relatively 
large, a reasonable estimate of the population variance can 


be obtained using the sample variance. The test is conducted 


with $ equal to .25, Ho equal to 75.25, and a equal to 


25) (9) 
lom 


a C* value equal to .0726 is obtained from Figure A-3 of 


or 1.875. Since the value of a is less than 3, 


Appendix A. Using this value of C*, the values of Cr. and 


Cp obtained from equation 4 are 


` 75.25 LI Оса 100075 -- 75.21 


Ω 
tt 


and 


tl 


С 75925 c 009735 = 75.20. 


R 


Since x is not between C, and C 


L R! the hypothesis that the mean 
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is between 75 and 75.5 cannot be accepted with a Type II 

error equal to .01. However, it may be possible to accept 

the hypothesis at a somewhat higher risk (8) of falsely 
accepting. Trying ß = .05, the new C* obtained from Figure 

A-3 of Appendix A is .35 and accordingly C, = 75.2 and C, = 75.3. 


L R 


Since x is between C, and C the hypothesis could be accepted 


L R' 
with a probability equal to .05 of falsely accepting. 
c. Example 3 

For the situation of Example 1, the decision 
maker wants to know for what values of a the hypothesis 
16.5 < u < 17 could be rejected. The answer can be found by 
estimating the maximum tail probability associated with the 
observed x using the procedure described in figure 2 with 


(х-н ууп 


а = 6.25 апа С* = рт рүе 7. Since a + C¥ > 3.5, 


the estimate of Pr would be Be Ф[.25- 7] = .227. Therefore, 
an a probability greater or equal to .23 would be required 


before the hypothesis could be rejected. 


B. CASE OF UNKNOWN POPULATION VARIANCE 
l. General Remarks 

A procedure for situations of unknown population 
variance which would parallel that of the previously dis- 
cussed variance known procedure cannot be developed without 
several stage sampling such as Stein's Procedure (reference 1). 
However, the previously developed methods can be extended 
directly to the case of unknown variance provided the null 


hypothesis is scaled by the standard deviation (see below) 
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a I 


and the noncentral t distribution is used in place of the 
normal. Details follow. 


a. A procedure was developed to test the hypothesis 


ΠΠ 


By < ro or equivalently < r where 


u is the true population mean, 


u is some specified value of interest near 
which the mean is to be tested, and 


r is the portion of the population o for which 
the test is to determine if the true mean 
is within ro units of Но: 

b. The procedure can be used to conduct tests about 
the displacement of the coefficient of variation (CV = ΠΠ 
from zero. 

2. Procedure Description 

The procedure for testing that — E rg Ss 
Erovided in figure 3. For testing that CV < су, figure 4 
should be used. Finally, the method for estimating Pn 15 
provided in figure 5. 

As the sample size gets large, the ratio = converges 
to the test value of r. To take advantage of eg 
the curves for a error level control are presented in terms 
of the value 7 instead of just C. This provides an indica- 
tion of the AS of r and the sample size n for which the 
с? unknown curves become close to those of the σ΄ known case. 


However, since the method of estimating P,, requires repeated 


T 
table entry and interpolation, the critical values presented 
in the tables of Appendix B are in terms of C directly. A 


Summary of the procedure follows. 
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a. For Testing That |u-u,| < ro 
Select the desired a or B error level. Consult 
the corresponding curves of Appendix A or tables of Appendix 
B. If the curves are used, multiply the value obtained for 


E- by the sample size (n). If the tables are used, determine 


n 
the corresponding value of C directly for the subject values 


of r and DF. Calculate the statistic 


> 


ا 
ρα.‏ 


If the test is to control the a error level then reject the 








hypothesis, with a Type I error equal to a, if the value of 
r is greater than the C value determined from the a curves 
or tables. For a test to control the 8 error level accept 
the hypothesis with a Type II error equal to B if the value 
of r is less than or equal to the C value obtained from the 
l- 8 tables. 
b. For Testing That CV < CV, 

To test the null hypothesis that the coefficient 
of variation does not exceed Ον. at the a level of significance, 
one need only apply a role reversal technique to the above 
procedure. The set equivalences 

(СУ < Су) = CS ον) = B^» we 

illustrate the fact that the hypothesis of interest is one 
of the alternate hypothesis treated previously, and can be 


managed by reversing the role and interpretation of a and 
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8 error level control from the previous test for 
- < š 
ΠΗ! E το 
1. For Controlling The a Error Level. Obtain 


the C value corresponding to the subject values of r (equal 


ل 
CV‏ 
О‏ 


the value 1-а. Calculate the statistic 


ο ) and DF from the table of Appendix B associated with 


^ 


r= 


хуп 


AZ 
σ 


If r < ο πεις ο Me nypochesits chat CV < cV with a Type I 








error equal to a. 
2, For Controlling The B Error Level. Obtain 
the C value corresponding to the test values of r and DF 
from the curves of Appendix A associated with the value 8 
(i.e., for R = .01, consult the curves for an a= .0l). 
Calculate the statistic г. If r » C, accept the hypothesis 
with a Type II error = 8. 
с. For Estimating P 


T 


The actual value of PT associated with the observed 


value of r will depend upon the true value of 
Se 


σ 


To obtain an estimate of the range of possible values for Pas 


the procedure detailed in figure 5 is used to calculate a 


value for Pr corresponding to r set at zero and r set at the 


test value (r). The actual value of Р will be between these 


two extremes. 
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Calculate 





For aerror control 


Go to the appropriate figure 

of Append. A for the associated 
alevel(Fig. 615) or tables 

of Append. B for parameter 

values not covered by curves 

of Append. A. Project vertically 
from the r scale to the appro- 
priate DF curve then horizcntally 
to the C/n scale. Multiply by 
үп Lo obtain C. 






If r > C, reject the hypothesis 
Ho with a probability = a of 


rejecting when in fact it is true 






For 8 error control 


Go to the appropriate table 
of Append. B associated with 
the desired 8 level(i.e. for 
8 =.1,go to the C value table 
associated with a Pr 1-8 =.9). 


Project down the corresponding 
r value column and across the 
respective DF row to determine 
the correct C value (Note: If 
the r value being tested is 
between two of the cable's 
columns ,interpolation is 
necessary.) 


If r < C accept the hypothesis 
Ho with a probability = 8 of 


: accepting when in fact it is 


false. 


Figure 3: Algorithm For Testing The Hypothesis и-н Š ro 





Calculate r= 





X Yn 





For controlling a error level For controlling 8 error level 


Go to the appropriate 8 table 
of Append. B (Tables 3 — 41) 
which corresponds to ang value 
equal to the desired a level 

ΠΠ e for as.] go to the: с 
critical value cable associated 
with 1-8 2.9). Project down 
the corresponding r value 
column and across the respec- 
tive DF row to determine the 
correct C value. [ Note : If 
the r value being tested is 
between two of the table's 
columns , interpolate for C by: 


πει lent с) 
High "Low 


+ 





If r < C, reject the hypothesis 
CV CV. with probability a 


of rejecting when in fact it is 
true. 


Figure 4: 
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Go to the appropriate a figure 
of Append. A (Fig.6-=15) which 
corresponds to an a value equal 
to the desired 8 value(i.e.for 
8 -.01, go to the a-.01] curves 
of figure 6 or 7). Project 
vertically from the r scale to 
the appropriate DF curve then 
horizontally to the C//n 
scale. Multiply by /n to get C 
[ Note : For parameter values 
not covered by the curves ,go 
to the tables of Append. B 
which provide C directly , 
instead of C//n.] 


Ifr>c , accept the nypo- 
thesis CV<CV with 


probability 8 of accepting 


when in fact it is false. | 


Algorithm For Testing The Hypothesis CV « CV 





Calculate r 





First calculate the estimate of P+} for r set equal to zero as follows: 





Go to Tables 3447.Pick a P-(i.e..5) and go to the associated table 


Project across the respective DF row and down the appropriate 
r column to obtain the correct C value. (Note: If the r being 
tested falls between two columns, ,interpolation 1s necessary. 









Go to a lower P+ value 
table(i.e. .2). 













as that table No 
been checked on a 


arevious step? 






Yes 


The Pr for this r 1s between the current table P+ value and the 


previous cable. For a closer approximation interpolate as follows: 


^ 


C 


: Гг 
ei) (Pr -P- )] 


low high “low high 





Now repeat the 
process to get 
an estimate of 
Pr for r = r. 





No 





estimates of P+ fo 
both r=0 and r=r bee 












The true P. lies between the two estimates for r=0 and rer 


Figure 5: Algorithm For Estimating Pp For o Unknown 
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3. Range of Parameter Values Considered 


Two general levels of r values, along with πες 
incrementation, were chosen to parallel the range considered 
applicable in the majority of situations. For example, in 
most cases where the test of interest is that Boon τσι 
r would range between zero and l. Conversely, most tests 
concerning the coefficient of variation would require that 
‚143 < Су, < 1 (or equivalently, 7 > r > 1). However, the 
critical values (C) for the total range of r (i.e., o to 7) 
provided in the curves and tables can be used for either 
test. The specific parameter ranges covered were: 

a. Curves 

Curves are provided for a error level control 
for both low and high values of r. 

l. For Low r Values (0 < r < 1). Curves are 
vided for a equal to .01, .025, .05, .1, and .2. 

ENG Values SS Curves are 
Provided for a equal to .01, .025, .05, .1, апа .2. 

b. Tabulated Data 

The same tables are used to present critical 

values for both a and ß error level control and for estima- 


ting P To obtain C values for ß error control, the tables 


т“ 
associated with a and PT equal to 1-8 are used. The values 


covered by the tables are a, P or 1-8 equal to .01, .025, 


má 
ar. DG O a 9 95, .975, anda .99. 
For each value of a, Επ, and 1-8, three tables are provided 


to cover the range of r from 0 to 7. 
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4. Example Applications 
a. To Test That Iu - uj SES 


It is necessary to determine if the output of a 
new system is between 11.75 and 12 units. The standard 
deviation in the output of the new system is believed (with 
confidence) to be between 2/3 and 3/4 of that of the old 
system (0 = 1.5). Using the conservative estimate of the 


improvement in the system deviation (i.e., 3/4), the test 


Te TS) 
"πο; - 


For the first 15 units, the values of the sample mean and 


parameter r would be set equal to dod Ges 
deviation were 11.68 and .3, respectively. The desired level 
of significance (a) is .05. 

The test 1s conducted with Ho equal to 11.875 and 
Е equal to 2.517. Using figure A-10 of Appendix A, n is 
determined to be .6, and accordingly, the critical US 
equals (.6) (¥715) or 2.324. Since r is greater than C, the 
hypothesis is rejected with a Type I error equal to .05. 

Dee Oo Test That Cv < CV. 

A relationship is known to exist between the mean 
and standard deviation of the purity of a certain chemical 
product. A new method of processing is being considered to 
improve the expected purity of the product. The decision maker 
15 interested in determining if the same degree of quality 
control would be obtained with the new process as with the 
current method which provide a standard deviation equal to 
.18 of the mean (i.e., a coefficient of variation equal to 


ELS). 
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Since a relationship between the mean and standard 
deviation is known to exist, and only 15 batches of the pro- 
duct have been produced by the new process, the population 
variance is considered unknown. The mean purity of the first 
15 batches was 4.6 parts per million with a sample variance 
equal to .25. The decision maker is concerned that if he 
does not convert to the new method, he will be rejecting, with 
a probability (a) equal to .05, an opportunity for a better 
product. The rejection region associated with a value of r 
181 to 5.56 (i.e., 9) LSeround by ¡consulting the critical 
value table with 1-8 equal .95. Interpolating between the 
table columns for r equal to 5 and 6, a value of C equal to 
17.69 is obtained. The value of r is calculated to be 35.63. 
Since r is greater than C, the hypothesis that the new process 
provides a quality control equal to or greater than the 
current method cannot be rejected with a Type I error equal 


po .05. 
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APPENDIX A 


CURVES FOR DETERMINING a OR 8 CRITICAL VALUES AND TAIL PROBABILITY 
(P 
T 





I. GENERAL 


The curves provided within were developed to simplify 
procedural calculations in lieu of the tabulated data of 
Appendix B. They do not cover the total parameter ranges of 
the data in Appendix B but were designed to maximize interpo- 
lation precision while covering what is felt to be the majority 
of applicable situations. For those cases which are subject 
to parameter values outside the range of the curves, Appendix 
B should be used. Figures A-1 through A-5 are relevant to 
the σ΄ known procedure. Figures A-6 through A-15 deal with 


the с unknown Situation. 
iL ACCURACY 


The accuracy achievable in using the curves vs. uSing the 
associated tables depends on the interpolation prowess of 
the individual user. However, to provide a general idea of 
the accuracy which may be expected, the following statements 


are considered germane. 


A. σ΄ UNKNOWN PROCEDURE CURVES 


a. C* Values For a or 8 Error Level Control 
Accurate interpolation to the third digit is possible 
m.e., C 


equals C* .005). 


* 
true interpolated B 
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lon Po Values 


Because of the semi log scale and interpolation 
between curves, the accuracy in the interpolated Р values 


cero, und + ‚005 гог 


is expected to be = + .02 for .1 < Р 


Р бИ [Гот Επι > .9. 


B. σ΄ UNKNOWN PROCEDURE CURVES 


a. Curves For The Range of r From 0 to l 
Interpolation accurate to the third digit is 
possible. 
b. Curves For The Range of r From l to 7 
Interpolation accurate to the second digit is 


possible. 
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C* Values For a Error Level Control W 
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C* Values For B Error Level Control With o Known; 
(2< a < 1.6) 


Figure Be 
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Tail Probability For o Known With 1.5 < a < 2.8 
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Figure A-8:  C//n Values For 
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Figure A-9:  C//n Values For o Unknown With a, Pq, or 1-8 = 
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Figure A-10:  C//n Values For c Unknown With a, P 
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Figure A-14:  C//n Values For o Unknown With a, P 
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C//n Values For o Unknown With α, 
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Figure A-15 
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TABLES OF a OR B CRITICAL VALUES AND P 





0,001 


3. 29056 
330076 
3,32634 
3,35758 
3,38453 
3,42115 
3,0616 
3,50459 
3,549281 
3,59654 
3,64434 
3,569228 
3З, 71193 
6.797 31 
3, 84091 
3, 89066 
3,9405 

3,9904 

u, 0031 
4,09031 
4.14928 
4.19027 
4, 24026 
4, 29026 
4, 34026 
4,39025 
4,44025 
4,49025 
4,54025 
4,59025 
4,64025 
4,63025 
4,74025 
ц, 79025 
4, 84025 
4, 89025 
4,94925 
4,99025 
5,04025 
5,09025 
5,59025 
€,09025 


Table B-1: 


0,01 


2,57583 
2,5085 

2,694 26 
2.525601. 
2,052 0% 
2,6834 
ο. τα 
2757382 
2,79867 
2,nu 223 
2 806251 
2,93409 
2,9r166 
34. 02995 
32, 07877 
ο ο πο 
3.17791 
3, 22701 
3, 27073 
ο 322663 
3.57653 
SHELA 
3.276142 
3,576 39 
3, 576 37 
3.62636 
3. 6 76 3f: 
ο, 52535 
3, 77635 
3.82635 
3.87635 
3.92635 
3497635 
4,02635 
4,07635 
u, 12635 
4,17635 
4.22635 
4,27635 
„32635 
4.82635 
5.32635 


0,025 


2,2414 
2,25254 
2,26629 
2.28519 
2,3093 
2,13711 
2. 36925 
2,40496 
2,44364 
2.40484 
uo eet 
2,57304 
2.61927 
2,666 52 
2.71453 
2 Feats 
2,81211 
2.86141 
2.91093 
2.9606 

3.01038 
3,0€024 
Зе ства 
3.16007 
2222008 
3,2601 
3.30999 
3.35990 
3,40997 
3,45927 
3, 503997 
3,559396 
3,602336 
3.65996 
3, 70996 
3.75996 
3, 80996 
3, 25996 
3.90996 
3.959396 
4,45996 
4,95996 


APPENDIX B 


0, 05 


1,95996 
1.959 52 
TU PIT 
1,997155 
2.01971 
2,04505 
2.07426 
2210699 
2,152865 
2.107108 
2, 22205 
2. 26549 
2.30996 
2,3928) 
2.40275 
2,45047 
2,4908 

2.45476 

259675 
2.64615 
2,695 73 
2.74544 
24 71952 
2,854511 
2.99502 
2,94436 
2.09492 
3,0449 

3,0948 
3,1447 
3. 1931166 
3, 2uu £ 
3. 294 86 
3. 364R6 
3, 39085 
3, цц? 5 
3,49425 
3, 54425 
3, 59425 
3, ЄЧЦЕЕ 
ц, 1408 5 
U,64UuRS 


0,075 


1.78046 


1,762335 
1,800237 
1, 21567 
1,83511 
1,P5851 
4, 8R565 
1,21625 
4950034 
1,9P666 
2,0258 

2,0671 

211020 
2,15491 
2,200R2 
2, 24772 
2,29542 
ашата 
2.33247 
2,U5158 
2,49094 
2,54049 
2,5901? 
2,63997 
2,689Р2 
2,73972 
2, 78965 
2,03961 
2,889 58 
2,233956 
2,098955 
3,039255 
3, 0P954 
3,13954 
3,12953 
23953 
3, 28953 
3,33953 
3,389353 
3,43953 
3,93953 
4,43953 


OF 


1,64485 
1.652306 
1.65327 
Toby ος 
1769667 
4127] 75 
Ἱ ντου 
er) 
1,893822? 
1.032875 
1,876 28 
4916142 
1,905735 
2.00147 
2.0464? 
2. 09252 
2139955 
5733 
2, 235€8 
ο, 2940 7 
2533353 
2, 30297 
203252 
2.487221 
2,5932 

d. SII ES 
2203175 
2.6P1F8 
2.134059 
2 70516 

2.83153 
2,2857 
2.331527 
2.98156 
3, 03156 
Qa 081 55 
313155 
3.121855 
3223155 
3З, 280165 
ο ο 
ц. 25155 


ο. 15 


1,43953 
1,44677 
1,45569 
1,1662 

4,882 22 
1,50367 
1,527F47 
1555214 
1,52159 
1,6136 

1.6003 

1.68557 
1.79855 
1,7661? 
1,8092 

1,852€8 
1, 899 33 
1,34592 
1,991522 
2,0415 

2,09708 
2243961 
2.10024 
2,2378 
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APPENDIX C 


PROCEDURE DERIVATION 
I. KNOWN POPULATION VARIANCE 


A. ERROR LEVEL CONTROL 
The a rejection region associated with the hypothesis 


that the mean is between dr and ар is defined by 


a s Pix < στ] + P[x > Cp] 


where the rejection region is chosen to be symmetric around 


uo (see figure C-1). Since the suprema will occur when 








-d=d 


1 £ È 
C=C, A uy dRo СЛН 


Ho 


Figure C-1: Density of x with u at ἃ, ( α -ΒΒ) 


u equals d, or ар, the critical region for the specified upper 


L 


bound on the Type I error (a) can be written as 


ОО = δες C. |given > ا‎ + P[x > C,|given ЕЕЕ e 
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Since the sample mean (x) will be normally distributed, 


this can be expressed in the form of equation (C-1). 


(C, -(u.- 8) ут (Cp-(u, 76) M 
0-1) α -ρ[χ-----ᾱ----]» [Ζ»----------] 
where z = 1 HN v N(0,1) 
Substituting Cr = Шс апа Ca = но +С, equation С-1 becomes 


(C-2) a spp z< ISCH J. pr z, LEN ı 


or 


a = ol SA - 0) + - ә 0 + 5] 


By reparameterizing, the number of parameters necessary 


for solving equation C-2 can be reduced from five to three. 


Letting a = an and C* = E equation C-2 becomes 
* * 
(0-3) a = 1+ ofa-C] - s[a4*€] 


το. το 


, equation C-4 is obtained and is used to define the 


Using the result that C 


τας 


= + 
Ho C, and 


Е = 





n 
rejection region. The hypothesis is then rejected, with a 


probability of falsely rejecting equal to a, if the sample 


mean is less than C, or greater than Св: 


L 
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I 
| 
| 


For each value of a there is a value of "a" such that 
the relationship between "a" and C* in equation C-8 is nearly 
linear (i.e., C* = K. +a) . For a > 1.7 the relationship is 
approximately linear (starting with 5 significant digits 
for the worst case with o = .25, and improving as "a" in- 
creases) with slope one for all values of o < .25. The value 
of Ko can be obtained by considering that if ó[a+C*] z 1, 
then equation C-3 simplifies to equation C-5. This can be 
rewritten to provide Kg = τι. When this linear relation- 


EDD occurs, Cr and Cn can be obtained using equation C-6. 


(C-5) a =1 + ofa - C - 1 or a= ofa - C 


B. P ERROR LEVEL CONTROL 
In order to determine an acceptance region for the 


hypothesis dr ST S «d with a probability no greater than 


р’ 
В of falsely accepting, the alternative hypothesis μμ] > ὁ 
1s tested. The a rejection region of the alternative hypothe- 
Sis will be the 8 acceptance region of the null hypothesis 

oe 


(1.е., dr a ἀπ) and will be defined by C 


L R 


(see figure C-2). 
The suprema will occur where u = dr or dac Thus, the 
upper bound on the Type II error (8) for the hypothesis 


d. < u < d, can be written as 


96 
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Figure C-2: Density of x with u at di ( 8 = shaped area) 


sup = Blu) = PIC 


<x< C„|given that u = u - ó] 
a 


1, 


This can be expressed in the form of equation C-7. 


(Cg-ug*6)fn (C -u tS) n 
G 


Re (Hr) = P[z < ] - Plze ر‎ 


[x -(u, -8) 15 


O 


where z = 


By substituting C, = Wn C and C, - TG, and noting 


L R 
that z will be normally distributed (i.e.N(0,1)), equation 


C-7 reduces to equations C-8. 


(c-8) oft.) = of (SHIH) әр (ΞΟΠ η 


πλ ο) ο. 


Q 
— 
EE 
— 
u 


Equation C-8 can be reformulated to resemble equation C-3 
which will permit use of the same solution technique as was 
used fora error control. Specifically, by multiplying 


equation C-8 by -1, adding 1 to both sides, and setting 
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l-a(H,) equal to a', equation C-8 reduces to equation C-3. 
The solution of this equation with a’ = 1-0 (H,) will also 
be the solution for a' = 1-6 (Hy) ° 

For each value of B there is a value of "a" such that 
the relationship between "a" and C* in equation C-8 is nearly 
linear (i.e., C* = Kota) - For a > 3 the relationship is 
approximately linear (starting with 3 significant digits 
for the worst case of B = .01, and improving as a increases) 
with slope one for all values of B » .01. The value for the 
constant (к) can be determined by noting that when 
Ф[а+С*] ~ 1, equation C-8 reduces to equation C-9. 
* 


1-8 = ola - oa or a-C = ү -8 ] 


* 


and C =a -K 


H 


(C-9) a 


zs 
E where К. = Ф [1 -8 ] 


The value of K, for 8 error control will be the negative of 


the value of Ky for a error control. To allow the same Ky 


~ 


value to be used for both a and 8 error control, equation C-10 
would be used for 8 error control when this linear relation- 
ship holds (i.e., a > 3 for this procedure). 


η (а- К.) Уп 


> Калы с een R = E 
(C-10) C. 7 Mg = απο те = 
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C. ESTIMATING TAIL PROBABILITY 


The actual value of p for an observed value of x will 


depend upon the true location of u within the interval dr 
to da (see figure C-3). Since u is unknown, the exact value 


Oy P,, cannot be calculated. By setting u = ug 7 8, the maximum 


possible value of P can be obtained from equation C-3 by 


“ΝΠ 
T for Ф, ара Еа АНЫ БАА 


T 


substituting P for C*. For large 


values of a+cC*, this reduces to equation C-ll. 








(C-11) P. = dla-c*] EOF атс > 3.5 
x Ho 
4 | ОР with Ημ. 
X d u i.” 
| L 0 UH) 
Рт with р = Ho” § = di 


Figure C-3: Dependence of P+ on Location of True u Value 


- 


39 
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II. UNKNOWN POPULATION VARIANCE 


A. TO TEST THAT |u- u,]| < го 
Ew Error Level Control 
The statistic 


O 


(x-u )/n 








72 
is used to perform this test and has the density of the 
absolute value of a non-central T, as depicted schematically 
in figure C-4. The associated a error level critical region 
15 defined when the non-central location is at r, resulting 
in a non-centrality parameter (X) equal to rn. The corres- 
ponding critical value (C) is obtained from equation C-12 


and the hypothesis is rejected, with a Type I error equal 


to a, if r is greater than C. 





0 A=rvn С 


Figure C-4: Density of the Non-central T(Absolute Value) 
Distributed f? 


(C-12) а = P[ P » C/ given л= гуп] 
2, B Error Level Control 


The suprema for the value of 8 will occur when 


A = rn. Equation C-13 is used to determine the value of C 
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which will define the acceptance region with a 8 probability 


of falsely accepting. 
(С=13) PCL | T | < C] = в 


By a slight reformulation, the same solution technique 
employed in solving for a critical values can be used here. 
Substituting a' = 1-8 into equation C-12 and solving for C 
will provide the critical value Season tony εν the desired 
level of a Type II error = B. This equivalence is shown in 


equation C-14. 
(C-14) в = РЕ IT| <c] = 1- PL |T| >C] or PLIT| >c] =1-8:=0 


Because of the equivalence between the C values for 


a and for 1-8, the critical values for a, P and 8 can be 


πι 
presented in the same table provided that for a desired 
value of a Type II error equal to 8, the table corresponding 


to a = 1-8 1s used. 


ШЕ LO TEST THAT CV < су 
1. 8 Error Level Control 
Tor -ormulate the procedure for testing CV < CV 
first consider testing the hypothesis E < r or equivalently 


CV 
HU 
| = 2| Ss ву Settlng Ho equal to zero, the null hypothesis 





can be rewritten as El < r. Solving equation C-12 with 


u = 0 andr = (хуп 


E | will provide the rejection region with a 





LOL 





λα. 


Type I error equal to a. However, rejecting the hypothesis 
that y < r, when r is greater than C, is equivalent to 
rejecting that CV > Ξ. This in turn is equivalent to 
accepting that CV < εν, with a Type II error equal to the 
value of a. As a result, it is possible to test the hypothe- 
Sis that the coefficient of variation is less than some speci- 
fied value (CV) by using the previous test procedure for 
testing that | < ro, provided the following differences 
are observed: 
a. Obtain the critical value (C) from the a tables 
(or curves) associated with a value of a equal to 
the desired value of a Type II error (8) and 


" 
with r equal to ον 


b. τετ is greater than or equal to this value of C, 
accept the hypothesis CV < CV with a Type II 
error = 6. 
2. a Error Level Control 
The same general formulation philosophy applies for 
Q error control of the test CV < СУ as for 5 error control. 


Again, first consider accepting the hypothesis = < r with a 


Type II error = B. This is equivalent to accepting that 


EV > = (or су). However, accepting the hypothesis CV >= 


with Type II error = 8 is equivalent to rejecting that 
i 
T 


Thus, the desired test that CV « CV, can be conducted for 


EV « (or CV) with a Type I error equal to the value 8. 


controlling the Type I error level by using the test procedure 
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for |u- uj] € rc and observing the following differences: 


a. 


B. ESTIMATING TAIL PROBABILITY (P 


Obtain the critical value (C) from the 8 tables 
(i.e., the table associated with the value 1-8), 
where the value of B is set equal to the desired 


value of a Type I error (a) and r is set equal 
dl 
CV, 


If P is less than this value of C, reject the 


to 


hypothesis that CV « CV with a Type I error 


equal to qo. 


т) 


The true value of P. will depend on the actual value of 


HAM 





to r (see figure C-5). The range of possible values for P 


ДЬ 


2S which, under the null hypothesis, may vary from zero 


Т 


can be estimated by calculating the Po associated with the 


observed value of r with A set at the two extremes zero and 


r/n in the solution of equation C-12. 






A=0 


Figure C-5: 


S 
S 


p(x at 0)- ⁄Z 


Pla at r/n)=|\\\\ 


2 





IT 
AT va. at wae 
κκ ЖАМАК 
A=rvn ` r 


Dependence of P. on the True Non-centrality Location 
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APPENDIX D 


SOLUTION TECHNIQUES AND COMPUTER PROGRAMS 
I. GENERAL 


The same general techniques were used in solving for the 
critical values in both the known and unknown variance 
Situations. This consisted of a combined proximity interval 
search and Newton-Raphson technique. A master program was 
used to monitor and control solutions for C values at various 
selected combinations of the respective parameter values. 
Solution accuracy was controlled by stopping the iterative 
procedure when the desired accuracy in the a or B error 
level was obtained. The specified accuracy for the case of 
known c? was .000001. In the case of unknown of, the 
specified accuracy was .00001 for all parameter values 


except the following: 


l. a= .01 with 125 degrees of freedom and r= .9 

2 a = .] with 90 degrees of freedom and r = .9 

3 a = .3 with 125 degrees of freedom and r = 1 

4. a= .3 with 60 degrees of freedom and r = 5 and 6 
5. a= .4 with 50 degrees of freedom and r = 2.5 


The specified accuracy for these latter parameter values 


was .00005. 
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ως σ΄ KNOWN SITUATION 


The programs used for this case were SIGMAKWN, CINCR, 
CINCRNEWRAPS, NEWTONRAPS, and NORMPROB. SIGMAKWN was the 
Master problem which stepped solutions along the selected 
а ог ß values for each specified "a" value. Controls for 
Switching solution subprograms were necessary for parameter 
values whose respective C* were too distant to permit solution 
solely by the Newton-Rapshon method. Equation D-l was used 


during proximity interval searching (CINCR and CINCRNEWRAPS). 
(D-1) ' F(a,a,C*) = 1- Ф[а * C*] + Ola = CF] 


Equation D-2 was employed when the Newton-Rapshon method was 


used. 


* * 
2.50663 (1-e[a*C. ,]te[a-C; ,]-a) 





* * F * 

(D-2) Es = С. m = ge + 
s 1-1 ж t ° 
oF / aC Exp[-.5(a+C, +) J+Exp[-.5(a-C, +) ] 


TII. с? UNKNOWN SITUATION 


The programs used for solving for critical values when 
g was unknown were MHUSIGRATIO, MHUSIGRATIOSOLENR, CINCR2, 
CINCRNEWRAP2, NCT, IB, and TRUNERR. MHUSIGRATIO was the 
master program stepping solutions along the selected r values 
for each specified DF. This permitted sole Newton-Raphson 


solution once the C for the first r value was determined. 
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To solve for critical values renee rated with testing the 
о? unknown two-sided composite null hypothesis, the equation 
P[r > C] =a (where r has a non-central T distribution) was 
solved for values of r between 0 and 1 using the method 
discussed in reference 2. This method requires use of 


equation D-3. 


2 2 
ον ο үа DF 1 
E 2 | (5-1) (n-1) 
(D-3) р[|т|>с] EA = ү y (1-у) dy] 
2 


where x =C + ( DF + С°) 


То solve the infinite series of the product of a Poisson and 
Beta distribution, the length to which the series must be 
summed in order to maintain a truncation error less than a 
Specified amount (c) was determined using the conservative 
estimate provided by equation D-4. Equation D-4 is solved 

by iterating m (the length of the series) until the left hand 


Side is less than e (this task is performed by TRUNERR). 


2. 
(D-4) [ “ pm Er 


mt 


This method of solving the non-central T is quite adequate 


for small degrees of freedom (DF < 30) and non-centrality 
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values (A < 3) but becomes lethargic for higher values. As 
a result, a second program was written in FORTRAN which used 


the nonecentral T subroutine of the IMSL packages. 
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[9] 
[101] 
[11] 
[12] 
[13] 
[15] 
[14] 
[16] 
[17] 
[187 
[19] 
[20] 
[21] 
[22] 
V 


V 
[1] 
[2] 
[31 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 
[112] 
[12] 
12] 
[15] 
[15] 


[1] 
[2] 
r3] 
[u] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 
y 


aGLOBAL INPUT;C(IPITIAÉ GUESS),DF 
JUNK+(SIGN+1),(LOOPF+O),(LL+30),(C+A-DFLC) ,(DELC+1 ) 
»L1xı((ALP20,05 )+(DF25))>1 

JUNK+(LL+60),(C+A-DELC) ,(DELC+10) 


L1:C+C+DELC 


+£5* 1 LL=LOOPF+LOOPF+i 
CI-LAM NCT C 
F+CI-ALP 

>L3 xı LOOPF z0 

-L3xiF»0 

SIGN- 1 


L3:-L1x10«SIGNxXF 


-LuxiDELC«0,51 


L2:C-C -DELC 


>[1 DELC«DELC+2 


LU YDERVIAT=-(CI-(LAM NCT(C-1E 8)))+1F 8 


+L2xıDEFRVIAT=O 
TLOXADELC«CPSDERVIAT) 
NEWTONRAPS2 

+0 


Poe "wexDIVERCEI?T AT *"3LL;'"LOOPS IN CINCENEWRAPI2(I="3I 3 ЛЕ у! 


C+A CINCR2 ALP 


NEWTONRAPS2;LL;CTEMP;JUNK 

aGLOBAL INWPUT:ERROR(DESIPED ACCURACY) 
JUI E-CLL-12),C(LOOP-0),C(DIVERG«*O) 
F+CI-ALP 


DI DEM BE 


C+-C-F+DERVIAT 

CI-LAM NCT C 

F+CI-ALP 

+0x1 ERROR>|F 

>52x1100P>1 

DERVIAT-(CI-(LAM NCT(C-1ER 8))):1Z72 8 


L2:>L1x1LL=L00P+L00P+1 


AER ENE ON MAPS DID NOT CONVERGES ITT JST; JL, SWITCHED 
DIVERG+1 

+Ox 1XILL #1 | 

C-CSTARUI;J-1] CINCR2 ALP 


ο ο πει στ ο οσο n REDE μα ο. 

AGLOBAL INPUT+¡PF,CYCLE(NUMBER OF LOOPS PFOUIRED TO 

a CONTROL TRUNCATION ERROR ) 

AGLOBAL FNPUT>P=DF+2,DF+-NS-1 

JUNK+( POS 70+(1+ж*1/м)), (TB2FEFRO*P га 0.5), (XX«0P4(D* ”+(0*2)) 
JUNK ( PPUPPER-POSSOX IBZFRO), (141) 

+OxX CYCLE =0 


Ll1:LL-POSSOxLÁ''*I 


PYTUPPER*PTUPPER+(P IB(I+0,5))x7.L 
POSSO+LL 
-—-Lixi(i1£e£CYCLE)zI-I41 
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V ARRAY-MHUSICBRBATIO;A;NS;,DF;P.C;.LAMV,CYCLE;,LAM;LA"V;LAMSQ2;XX; LAM 
[11 aCLOBAL VARIABLES=LITTR(VECTOR TEST inpusstate LITT) 
[2] aDFVECT(VECTOR OF DEGREE OF FREZDOH);ALP- LOS; EPROR=DESIRED PRECISION 
[3] JUNK-(I-1),(A-LITTRUÜ11),(KXILL-0) 

[4] DFVECT-,DFVECT 

[5] CSTAR-((1*pLTTTR )x (oDFVECT))40 

[61 CSTAR+((pDFVECT),(1tpLITTR))pCSTAR 

[7] PRECISION-A*ERROR 

[8] L1:DF-DFVECTÜI] 

[9] "РАТА FOR DFPzs!;DF;' WITH ALPHA=';ALP 

luc] —* LOOP TIME LITTR . ο 

[11] P+*DF+2 

(121 LAMV-LITTRx(DP*1)*0,5 

[13] CYCVECT+LAMV TRUNERR ERROR 

Γι] CYCLE+CYCVECT[A] 

[15] TIMEO+121 

[16] LAM-LAMSQ2(1] 

[17] >L5x14>0 

[18] А=0,0001 

[191 L5:C-A CINCRNEWRAP2 ALP 

[20] +L6x \DIVERG=0 

[21] C+A CINCR2 ALP 

[221 Γδις2 

[23] CSTAB[I;1]-DFVECT[I] 

[2] L11:-D513xXCSTAR(I;Jls 1 

[25] +L7xyJ=2 

[261] TI"FO-rI21 

[27] CY^7LZE-CYCVECT(J -1] 

[28] LAM+LAMSO2( 7-11] 

[29] CI+LAY ПСТ С 

[30] DERVIAT-(CI-(LAM NCT(C-i1E &)))#1" 8 

[31] NEFTONRAPS2 

[321 >LE7xvDIVERG=0 

[33] KILL-1 

[34] C-CSTARÜI;J-1] CIPCRNEWRAP2 ALP 

[35] KILL+O 

[36] L7:CSTARU(I;J]l-C 

[37] TIMEJ+*((121)-TIMEO)+60 

35] '  LOOP,' 'es( TIMES; ! F'eLITTREJ-1];' eC 
[39] >£13x1J=2 | 
[40] MULT-1 

ul] L15xii1s((LITTP[J-1]50,50 )4((J-1)spLTTTP)) 

[42] MULT-(LITTR[J]-LITTRI[J-11)s(LITTRU[J-1]-LTTTPUJ -2]) 
C43] L15:C-C«"ULTx(C-CSTAR[I;J-11) 

[uu] L£13:+£11x1(2+pLTTTR)xJeJ+1 
ν.Ι..'.-------------------------- еске еы A ш кишене ' 
(467. >£1x3(1+pDFVECT)x=I*I4+1 

[47] TEMPLITTReO,LITTR 

Για] APRAY+,CSTAR 

[09] APRAY*TEMPLITTR, ARRAY 

[50] APRAY+*(I,(7-1))pAPRAY 


tIl 





[9] 

110] 
[11] 
[12] 
[13] 
114] 
115] 
[16] 
[17] 
fis) 
£19] 
[20] 
[21] 
[22] 
23] 
[24] 
[25 ] 
[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[321 
[33] 
[34] 
835] 
[36] 
[37] 
[38] 
[39] 
C40] 
C41] 
Γι 21 
C43] 
[+] 
[+5] 
C46] 
[47] 
[u P] 
[49] 
[501 
[511 
[52] 
[531 
[54] 
[55] 
[56] 


V APPAY+MHUSIGRATIOSOLENR OLDCAPR SASNS 397 3PsCs3CYCh? 3 LAMV 3 LAMSON SYN LLY 
ΠΡ AR GAL LES =L0T TRU ECTOO απο умт /STEis LIT) 
ADFVECT(VECTOR OF DEGREE OF FRPEDO");¡ARPSTIOS:FRROP=DFETPED PPECISION 
AOLDCAP2(2xpDFVEC™ APRAY OF PREVIOUS C*S YTTE ROVA POVAT, τη 
GOMOD IS AND 1 CO 1 EQUAL TO O,DFVFCT) 
ACESTIMATE( VECTOR OF ESTIMATED Ç POP FTRET ERTMENT OF LI™?TA IM EACH DF 
CESTIMATE*LITTR(1] CESTIMATOR OLDCARR 
DPVECT-,DFVECT 
JUNK-(I«-1),(A-CESTIMATEC1]1),C(KXTLL-0) 
CSTARe((1+pLITTR)x(pDFVECTINO 
CSTAPe-((oDPVECT).(1*nLITT?))pCS7AR 
PRFCISION+1+ERROR 
L1 :DF+DFVECT[I] 

'DA7A FOP DFzs';DP;' WITH AKPHA=!;ALP 
' LOOP TIME LITTR С! 

P+DF #2 
LAMVeLITTAX(DF+1)*r0,5 

CYCVFRCTELANV TRUNERR ERROR 
CYCLE-CYCVECT[1] 

TIMEO+121 

LAM-LAMSQ2[1] 

CeCFSTIMATEL IT) 

CI-LAM NCT C 

DERVIAS-(CI-(LAM NCT(C-12/8)))*1E ^8 
NEWTONRAPS2 
L6 :Je2l 

CSTAR[I;1]-OLDCAPR(I*1;3] 
D11:-—D513*4CSTAPU[I;J1s 1 

+L7xı1J=2 

TIMEO-I21 

CYCLE-CYCVECT[J -1] 

LAM-LA"SQ2UJ -1] 

CI-LAM рст C 

DERVIAT-(CI-(LAM NCT(C-1E 8)))saiF 8 
NEWTOFRAPS2 

>L7X1DIVERG=0 

KILL-1 

C-CSTAR[I;J-1] CINCRNEVRAP2 ALP 
KILL+0 
L7:CSTAPR[I;J]-C 
TIMEJ+((121)-TIMEO)+60 


' лел еы TIES LITT ету! p 
MULT+1 

>»L14xıJ#2 
MULT-(LITTR[Z]-LITTRUJ-1]):(LITTR[J-1]-OLDCARRI1;3]) 
>L15 


Li15:-Li15xiis((LITTP[ZJ-1]30,8)«-€((7-1)20LTTTR)) | 
MULT-(LIP"TR(J]-LITTR[ZJ-1]):(LIT7TR[J-1]-ZTTTPRUJ-2]) 
L15:C-Ce"ULTXx(C-CSTAR[I;ZJ-1]) 
L13:2L11xi(24pLITT"?)zJ-J*41 

CSTAR[I;1]-DFVECT[I] 


-Lixi(1topDFVECT)2I-I-*1 
NREMPLINTREN,LITTR 

Minn YG ας 

ARRA Lehr hE OR ARRAY 
ARRAY-(I,(0-1))0APRAY 
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V M-LAMV TRUNERR FRROP; JUNK SELAM:N:K:ERR 
t1] AGLOBAL INPUT: LITTR,ERROR 
[2] TUNK+-(K+1), (4-10) ,(LAMSQ2+0,5*XLAMVx«2) 
[3] L1:-L3xiLITTR(EÉE]»1 
Γι] -Li11xiLAMVIK]s6 
(57 ΓΙ (11), (Ο5ΠΑΡΓΙ Κι ]ς 1) 
[6] L11:-L3x1LA"VVCX]20,1 
[7] ier 
(8] >L4 
[91 L3:NH-((CFZAMYVCX])*2) *2x[ LA AMYEX1)- 2 
[101 FLAHM-LAMSQ2[K]x2,71828183 
[11] L2:ERR-(ELAMtN*T1)*N«*1 
[12] NeN+1 
[13] >L£2x1ERR>ERROR 
(145) Lu:M-M,97-1 
(15] >£1x1(1+pLAMV)=K+*X+1 
. V 


Q REB IB Q:VsMA5SXS5RSBR.H:QQ I ;P.J 5 РРЕСІБ 
[11 AGLOBAL INPUT; XXCDFA(DF*C*2)), PPFCISION( EQUALS 1+ERROR) 
[2] X-V/WXxV-XXs0,5 
[31 J+(T+0),(R+0) 
C4] L1:4-(050 0) «(0-LO)xQzl Q 
E I+I+1 
[6] >L4x1(+/V)=0 
(7], M-[PTI(aGaPRECISIOF*(P!A*P-1)x(T/X)xP)*0,6931571806 
[81 В+0 
[91 >[3х\А=1 
[101 BBe(i1M)o,24M 
[11] B«-(M,M)pCCuyM) -A) xi M 
[12] B+x/(BxBR)+1-BB 
[13] Be-(H-(XxP)xP!A*P-1)9 , X RXP& Pt AM 
[14] B+H+r/BxXo,xıM 
[15] £3:f+(OQ, Q0 pC Q+1-100 )tP+0-190+L0 
[16] BB-(1900)25,210Q 
[17] BRex/(HxBR)+1-BB 
(18] BB<(((Q!'P+0-1)xX*P)e,xBR)x(1-7X)°,*Q- OQ 
[19] B+B++/BB 
[20] R+R+(V\BxI=1)+(V\1-B)xI=2 
[21] Lu :F#<P 
[22] J<(V<<—V),(O<H),( P<0 ) 
[23] Xev/1-XX 
[241 >£51x1(0<pX)aAlT=1 
y 


v CESTIMATE-NFURA CESTIMATOR PREVIOUSAPR.LASTO.JUNK;I;DFOLD;MULT 
[1] aGLOBAF INPUT-DFVECT ASSOCIATED WITH NEW LITPR 
[2] TUNK+<(R+<1 5p PPEVIOUSARR) ,CI-1),CCESTIMATE-10) 
[3] LAST2*(R, 2)*PREVIOUSARR 
[u] DFOLD<, 1 +PREVIOUSARRT ; 11 
[ 5] MUL Pe( NEWRY HAS TILL: 21) +(LAST2[1; 2]-LASTO[ $511) 
[6] L1:>L2xı1=(+/(IFVECT=SDFOLDET))) 
[7] СЕ$ТТМЛТЕ&+СЕ5$ТТМАТЕ, (Ми Тх(145Тт2[1+\1:2]-1А5Т72[1+1;1]))+145тТт2{Т+1 27 
[8] Z72:>+Lixı(1+pDFOLD)zI+I+1 
V 
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SYNTAX: PTAIL+DFVECT NCT CVECT 


PERIS FUNCTION CALCULATES THE UPPER AFD LOWER TAIL PROBABILITY 
OF A NON-CENTRAL T DISTRIBUTED R,V, RELATIVE TO A DPSIRFD 
VALUE C(FOR EXAMPLE,PTAIL zPROBE|R,V, I»C1),DFVECT IS A VECTOR 
SCALAR) CONTAINING DEGREES MF FREEDOM OF INTEREST,CVYECT IS A 
VECTOR(OR SCALAR) CONTAINING C VALUES OF INTEREST,FUNCTIONS ASSOC- 
ο ο WITH(AND CALLED AUTOMATICALLY) APE TRC INCOMPLETE RETA) AND 
TRUNERRCDETERMINES LEVEL TO WHICH THE INFIFITE SEPIFS ASSOCIATED 
ο. τς NCT CDF MUST BẸ SUMMED IN ARDER TO RFSTRICT TRE TRUNCATION 
DEOR TO THE SPECIFIED VALUE).THFE CPU TIME REQUIRED TO EXECUTE 
ΠΩ FOR ONE VALUE IS GFEFATLY FFFFCTED BY TÜ"E SPECIFIED 
VALUES OF DF,ERROPR,C,AND ACAS DEFINED BELOW),TRF LARGER THE DF AND 
A,AND SMALLER THE ERROR THE GREATED THE CPU TIYE,AS A POUGE IDEA- 
BORZERROR+,0001,.01sAs,2,1sDFs16, ,O1sCPUS2 SEC, PEP VALUE,FOR 
AE 000001,.284S1,15SDFS50,2<CPUS120 SEC, PER VALUE,TO COPY 
ANDY IHR FUNCTIONS NECESSARY FOR FCT COPY GROUP FONCENTRALT 


GLOBAL LUPUT | 
1,ERROR+DESIRED LEVEL OF ACCURACY(T.F, FPROR<,0001 WILL PROVIDF 
RESULTS WITH A TRUNCATION ERROR OF LESS THEN ,0001 FROM THF 
TRUE TAIL AREA) 
2 A+LOCATIOF(DISTANCE) OF TRE ИСТ FROM ZLRO(NOTE:THTS DETERMINFS. 
THE NONCENTRALITY PARAMETER -LAM-), 
3,0PTION<O0 OR 1(0 WILL READ IN GLOBAL -A- APD CALCULATE -LAF- INTERNALLY 
„1 WILL READ IN -LA"- GLOBALLY AND CALCULATE -A- INTERPALLY FOR 
EACH DEGREE OP FREEDOM 
4 LAMeNONCENTRALITY PARAMETER=AxX(DF+1)*,S[NEED ONLY BE SPECIFIED IF 
OPTION+1] 
S HEANING*O OR 1(0 WILL SUSPRESS PRIFTING OF A DESCRIPTIVE READING 
ABOVE THE OUTPUT PTATITL ARRAY) NOTE: SPACING FOR TEE DFSCRIPTIVE 
HEADER IS DESIGNED FOR )DIGITS 4,5, OR 6, 


OUTPUT 
OUTPUT FROM THE FUNCTION WILL CONSIST OF AW APRAYCSCALAR) OF VALUFS 

wie CVFC!L AS THE FIRST ROW, ,DFVECT AS THE FIRST COLU*"N.-A- AS THE 

SECOND COLUMN,-LAM- AS THE THIPD COLUMN, AND THF AS SOCIATED PTAIL IN 


THE RE'"AINDING COLUMNS AND ROWS 


V CYCLE-LAM TRUNEPR ERROR;JUEK;ELAM;N;ERR 
[1] JUNY-(CYCLE-CCL LA) 2) e3xLLAV) , C LAV802-0, 5x LAM x2) 
12] FLAM+LAMSO2%2,71828183 
[3] L1:EPR-(ELAMSCYCLE*1)xCYCLES*1 
[u] CYCLE+CYCLE+1 
[5] >L1X1ERR>ERROR 
[6] CYCLE-CYCLE-1 

g 
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Y PTAIL-DEVECT NCT CVECT,JUFK;DP;C,P LAMESA? XX; LL; N; IBZERO; PT 
[1] JUNK+(I+1),(DFVECT+,DFVECT) ,(CVFCT+,CVECT), ( PPECISION+1+%ERROR ) 
[2] PTAIL-0,0,0,CVECT 
[3] L3i1:J7UNK-(P-0,5xDF),(DF-DFVECT(I]) 
am -L2xi0PTION-1 
[5] LAM-Ax(DP*1)*0,5 
[61] -L3 
[7] L2:(A4-LA"*(DF*1) 
[8] 213:PTAIL+PTAIL,DF,A,LAM 
[91 CYCLE-LAM TRUNERR ERROR 
[101 Lu: JUNK+=-(PO0SSO+1++LA"SQ2),(J+1) 
[11] L5:JUNK-(IBZERO«-P IB 0,5),(XX+DFs(DF+CVECTLJ]#2)),(N+1) 
[12] „о mno, (LL+POSSO) 
[13] >L7xıCYCLE=0 
[15] L6:LL-LLxLAMSQ2*N 
[15] PT7-PT«(P IB(Nt0,5))xLL 
[16] >26x1(1+CYCLE)*NeN+1 
[17] L57:PTAIL-PTAIL,PT 
[18] >25x1(1+pCVECT)*J+*J+1 
[19] >251x1(1+pDFVECT)=I<I+1 
[20] PTAIL-((1*0DFVECT),(340CVEC"))p PTAIL 
[21] >0x1HEADING=0 | 
moo] ° DF | А | LAM |+--C AND ASSOCIATED PTAIL --+! 


V R-P IB OsVishsA3sX3:B3PRB3H3003T3P3sJ3PRECIS 
[1] a INCOMPLETE BETÀ FNC,- REQUIRES DATA VECTOR XX 
[2] X-V/XXxV-XXs0,5 
πο] J+(I+0),(R+0) 
[tu] L1:4+(Q0=10)+(0Q- л 
R5] I+I+1 
[6] »Luxı(+/V)=0 
[71 M+<4[ P[ (@ePnECTSION:(P'A+P-1)x([ /X)*P):0,653931571806F 
[8] R+Q 
r9] +L3xı14=1 
[10] ZR+(1Y)o,>2ıM 
[111 А+(м, 7 )р( (М) -Ад) +1 
[12] B+x/(BxRR)+1-BR 
[13] B+(H*-(X+xP)xP!A+P-1)o,XRxP+P+ı 
[15] BeH++/BxXo, *. 1 
1157.23:Re(00,00)0(0+1-100)*+P+0-100+1 0 
[16] BB+*(100)0,2100 
[17] BBex/(HxBR)+1-BB 
[181] BR+(((Q!P+0-1)xX*P)o,xBB)x(1-X)0,*r0-100 
[19] BeB++/BR 
[20] R+R+(V\BxI=1)+(V\1-B)xI=2 
[21] L5:H«-P 
[22] J«-(Ve-V),CO-H),CP-Q) 
[23] X+V/1-XX 
[24] -Li1x1i(0O«pX)AI-1 
V 
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